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We investigate the bulk orbital angular momentum in a two-dimensional time-reversal broken topological 
superconductor with the Rashba spin-orbit interaction, the Zeeman interaction, and the 5-wave pairing potential. 
Prior to the topological phase transition, we find the crossover from 5 wave to p wave. For the large spin-orbit 
interaction, even in the topological phase, L z /N does not reach -1/2, which is the intrinsic value in chiral p- 
wave superconductors. Here L z and N are the bulk orbital angular momentum and the total number of electrons 
at zero temperature, respectively. Finally, we discuss the effects of nonmagnetic impurities. 

PACS numbers: 74.20.-z,74.20.Rp 


I. INTRODUCTION 

Topological superconductors (SCs) have been intensively 
studied due to their fundamental interest in the physical real¬ 
ization of Majorana fermions^ and their future application to 
topological quantum computation 2 . A typical design princi¬ 
ple for two-dimensional time-reversal broken topological SCs 
is to realize a spinless chiral p -wave SC 3 . Possible candi¬ 
dates are the surface of a topological insulatoi® and a quantum 
anomalous Hall insulator in proximity to 5-wave SCs 5 . More 
promising candidates are conventional spin-orbit coupled sys¬ 
tems with the Zeeman interaction and the 5-wave pairing po¬ 
tential®. In fact, their one-dimensional analog was theoreti¬ 
cally proposed 10 11 and later experimentally realized in semi¬ 
conductor nanowires contacted with v-wave scPEi. 

In a chiral /7-wave SC, each Cooper pair carries the orbital 
angular momentum (AM) l - 1, which leads to the Chern 
number C = 1 and the bulk orbital AM L z /N = 1/2 at zero 
temperature 16 25 . Here N is the total number of electrons, and 
we set ft = 1. Note that L z /N is known to deviate from 1 /2 due 
to retroreflection by edge roughness 23 and the opposite chiral¬ 
ity with the vorticity w_i = +2 induced around the edge 2627 . 
However, these effects do not occur without the edge, disor¬ 
der, or a vortex, and it is natural to call L z /N = 1/2 the in¬ 
trinsic value in chiral /7-wave SCs. In higher-order chiral SCs 
with l > 1, the Chern number and the bulk orbital AM are 
equal to C - l and L z /N = £/2, respectively. Although 
it was recently found that the bulk orbital AM is totally sup¬ 
pressed due to the edge states without the particle-hole sym¬ 
metry by themselveJS®, still we can call L z /N = 1/2 the 
intrinsic value in chiral SCs. Thus, the Chern number and the 
bulk orbital AM appear not to be independent of each other. 

Generally speaking, topological invariants are robust 
against symmetry-preserving perturbations which do not close 
the gap, while the orbital AM is not conserved in the 
presence of spin-orbit interactions (SOIs). In the specific 
model with the Rashba SOI a, the Zeeman interaction h, 
and the 5-wave pairing potential ^P®, which is dubbed the 
Rashba+Zeeman+v-wave model below, the Rashba SOI nei¬ 
ther closes the gap nor changes the Chern number. One ques¬ 
tion is whether the Chern number and the bulk orbital AM are 
not independent even in the presence of SOIs. If not, a sub¬ 


sequent question is what the bulk orbital AM characterizes in 
topological SCs. 

One possible answer to the second question is the impurity 
effects. In contrast to the Chern number which characterizes 
a topological SC, the bulk orbital AM may characterize how 
close it is to a spinless chiral SC with the same Chern number 
and how fragile it is to nonmagnetic impurities. Topologi¬ 
cal SCs are not always robust against nonmagnetic impurities, 
namely, the midgap bound states appear and the critical tem¬ 
perature T c is suppressed 31 33 . These impurity effects in the 
Rashba+Zeeman+v-wave model depend on h and are similar 
to those in a chiral /7-wave SC for large h. 

There is another evidence that the bulk orbital AM is related 
to the robustness against nonmagnetic impurities. The impu¬ 
rity effects on the bound state in a single vortex core were 
investigated in the Rashba+Zeeman+v-wave mode® There 
are two types of vortices: a parallel (antiparallel) vortex when 
the signs of the vorticity and the Zeeman interaction are op¬ 
posite (same). Note that our sign convention of the Zeeman 
interaction is opposite to that in Ref. [34| The authors found 
that the low-energy scattering rate is suppressed in the case of 
an antiparallel vortex compared to the case of a parallel vortex, 
and that the robustness of the bound state against nonmagnetic 
impurities depends on the Rashba SOI. Thus, the impurity ef¬ 
fects may be characterized by the bulk orbital AM rather than 
the Chern number. 

In this paper, we investigate the bulk orbital AM in a two- 
dimensional time-reversal broken topological SC. We focus 
on the Rashba+Zeeman+v-wave modef®. We calculate the 
bulk orbital AM both by the Berry-phase formula 25 ^ and in 
the circular disk. One advantage of the Berry-phase formula 
is that the wave functions in the reciprocal space are simpler 
compared to those in the real space. We explain the effects of 
nonmagnetic impurities in the absence 31 33 and presences of 
a vortex. 


II. MODEL AND METHOD 

The Rashba+Zeeman+v-wave model 6 ® is represented by 
1 r d 2 k t 

H ~2 J ^ ^ kT7 + — k y cr x )T z 
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(1) AM L z are calculated by 


where TV = ^ ,-c\ ] T is the Nambu spinor, and 

K k\ n -ki ~k\ 

cr’s and r’s are the Pauli matrices for the spin and Nambu 
spaces, respectively. & = k 2 /2m - p is the kinetic term in 
reference to the chemical potential p. There are four dis¬ 
persions, ±E k + = ± ^jg 2 + (ak) 2 + h 2 + A 2 + lg 2 with g 2 = 

yjg 2 (ak) 2 + % 2 h 2 + h 2 A 2 . Therefore, the gap between in¬ 
closes at k = 0 for h c = ^p 2 + A 2 , above which this model 
has the nontrivial Chern number C = -l 7 9 . It can be eas¬ 
ily understood by representing Eq. (JTJ) in the band basis 9 . In 
Fig.pTa) for h < h c , the gaps open due to the interband s-wave 
pairing potential A K k = hA/f k , the intraband (p-ip )-wave one 

A(f-^) = iake~ l(f> A/ f k for the lower band E^- — %k ~ fk* and 

k 

the intraband (p + ip )-wave one A^ +ip ^ = iake l(f) A/ f k for the 

k 

upper band E ok+ = £ k + f k , where f k = ^/(ak) 2 + h 2 . For 
h> h c , the upper band is far away from the chemical potential 
as shown in Fig.[ljb), and the model is mapped onto a ( p-ip )- 
wave SC with C - - 1. Below we set m - 0.5, p = 0.5, and 
A = 0.35, leading to h c = 0.61, and change a or h. 

First, we calculate the bulk orbital AM in the reciprocal 
space. The wave functions with the dispersions +E k -, +E k+ , 
- E k+ , and -E k ~ are simply given by 


c = E ( > - 7)[(*dy - viy) - (u 2 0i j - V 2 0ij )l (3a) 
ij 



Lz= ~f Z (/ ■ /)( '4- ■ v *.-A ( 3d > 


respectively. Note that the total number of electrons A as well 
as the bulk orbital AM L z changes by changing a or h with p 
being fixed. Nonetheless, the bulk orbital AM per electron 
L z /N plays an important role when we discuss the impurity 
effects below. 

The bulk orbital AM can be calculated not only in the re¬ 
ciprocal space but in the real space. We consider the circular 
disk with the radius r c . The Bogoliubov-de Gennes (BdG) 
equations in the real space are given as 


[H 0 (-N)t z + At x ] 


u(r) 

v(r) 


uir) 
Hr) 


(4) 


where 


xh 2 - g 2 k -hE k -)(e k -g 2 k +tkE k -)- 
-iake^ih - &)(£* - g\ + ZkE k -) 
-A(h - { k )(h 2 - g 2 k - hE k J) 
-iake" 1 ’ A(h 2 - % 2 ) 

-iake-^Qi + &)(£? + gl + frEk+Y 
(h 2 + g k + hE k+ )(£ k + g k + XkEk+) 
iake- ,,p A(h 2 - g 2 ) 

A (h + % k )(h 2 + gf + hE k+ ) 


N k+ 


N k - 


-A(h + £k)(h 2 + g 2 k + hE k+ ) 
iake itp A(h 2 - £?) 

(h 2 + g 2 + hE k+ )(£ 2 + g k + £ k Ek+) 

. iake ilp (h + &)(£? +g 2 k + £ k E k+ ) . 

—iake~" ! 'A(h 2 - £ 2 ) 

A(h-{ k )(h 2 -g 2 k -hE k .) 
iake-'Hh-mtl-gl+tkEk-) 
Xh 2 -g 2 k -hE k -)(e k -g 2 k +tkE k -). 


Ukll 
ie l<p u k 2 \ 
~Vk22 
-ie^Vk i2. 


(2a) 


-ie-'^Uku 

u k22 

ie~ i,p v k 2 1 ’ 
-Vkll . 

(2b) 


v k u 

ie'^vn i 
Uk22 

ie ltp u k {2 


(2c) 


-ie^Vku 

Vk22 

-ie-^uta i ’ 

UkU 

(2d) 


respectively. Here N k + are the normalization constants, and 
u k ij and Vkij are real functions of k only. In the polar coordi¬ 
nate, the Berry curvature for the wave function | ug n ) is defined 
by = i((d k U£ n \d(pU£ n ) - (d(pug n \d k ug n ))/k, and the integrand 

for the bulk orbital AM is given by (A^ x k) z = -i(u^ n \d^u^ n ) 9 
where A^ = i{u^ n \d k .u^ n ) is the Berry connection 25 . Since the 
Berry connection can be interpreted as the expectation value 
of the position operator x in the reciprocal space, the integrand 
is interpreted as that of x x p. The Chern number C, the total 
number of electrons A, the total spin S z , and the bulk orbital 


tfoHV) 


£(r 9 6)-h —aL-(r, 6) 
aL+(r, 6) £(r,6>) + /* ’ 


(5) 


with f(r,0) = -(d 2 /dr 2 + (l/r)d/dr + (l/r 2 )0 2 /d6 2 )/2m - p 
and L+ = e ±l6 (d/dr±(i/r)d/d6). With the use of the rotational 
symmetry along the z axis in the circular disk, the solutions 
can be expressed as 


u n (r) 

- p ine 

« T (r) 

e w ui(r) 

Vn(r) 

— c 

V| (r) 



-e' fl i’|(r) 


Here n is the quantum number. When [u n (r) T , v n (r) T ] 
is a solution with the energy E n , [uZ n _ 1 (r) T , (^*) T ] = 

[~icr y v n (r) T ,icr y u n (r) T ] is the solution with the energy -E n . 
We use the normalization condition, 



rdrd#(|w 0 -(r)| 2 + \ V(r {r)\ 2 ) = 1, 


(7) 


and the boundary conditions, 


a 

d~r 


Wf(r) 1 
(r) 
vj, (r) 

_v r( r )J r=0 


= 0 , 


w T (r c )' 
U[(r c ) 
vj.(r c ) 
-vt(r c ). 


= 0 . 


(8a) 


(8b) 
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FIG. 1. Effective pairing potentials induced by the 5-wave pairing potential. Solid and broken lines indicate the electron and hole bands, 
respectively, (a) For h - 0.33 < h c , the pairing potentials A ( _f +ip) open the gaps in the lower and upper bands, respectively, (b) For h - 0.67 > h c , 

k " 

the upper band is away from the chemical potential, and A { ^~ ip) opens the gap in the lower band. We set a = 1. 


The bulk orbital AM in the real space is expressed as 

n c /-*r c 

L z =2n V V I rdr[(n\u\(r)\ 2 

ilTo J o 

+ (n + \)\u{ (r)| 2 )/(£f) - ((« + l)|v'(r)| 2 
+ n\v\ (r)\ 2 )f(-E?)l (9) 

with the /-th eigenvalue E with fixed n and the Fermi-Dirac 
distribution function f(aj) = + 1). Note that we focus 

on f = 0. In order to solve the radial BdG equations, the 
second-order finite difference method with N g real-space grid 
points is used. Thus, the BdG equations become the 2N g x 2N g 
matrix eigenvalue equations. We set the number of the grid 
points N g = 1000, the cutoff n c = 512 - 1, and the disk radius 
r c = 100. The unit of the real space is defined by the lattice 
spacing of the tight-binding model in Ref. [33] 


III. RESULTS 

Figures [2] and [3] show the h and a dependences of C, S z /N, 
and L z /N. L z /N obtained by the Berry-phase formula com¬ 
pletely coincides with that in the circular disk. For small a , 
we obtain L z /N ^ 0 for h < A and L z /N - -1/2 for h > h c 
as shown in Fig. [^ a ). These are expected from the v- and 
(p - /p)-wave behaviors, respectively. Remarkably, in the in¬ 
termediate region A < h < h c , we obtain the nonzero L z /N. In 
contrast to C, L z /N does not jump but shows the crossover. As 
a increases, the region with L z /N ^ 0 gets narrower. For large 
a , as shown in Figs.[2jc) and[2jd), L z /N never goes to -1/2 
even in the topological phase. In Fig. [3] we find that a tends to 
suppress L z /N. Although the Rashba+Zeeman+v-wave model 
is mapped onto a chiral p-wave SC, these two are different in 
terms of the bulk orbital AM. These results are summarized in 
Tab. [D 

Let us discuss the origins of the h and a dependences of 
L z /N. In the band basis, the interband s- wave pairing poten¬ 
tial A^ opens the gap at k = J2m/j and Eok+ = —Eok- = 
a/2 ma 2 ji + h 2 as shown in Fig. Ill Therefore, v wave is dom¬ 
inant for Eok+ < A, namely, h < h* = A 2 - 2ma 2 jj. In the 


TABLE I. Summary of the h and a dependences of the bulk orbital 
AM per electron L z /N. For small a , h* = A. 


Small a 

Large a 

Chern number C 

h < h* 0 

Finite 

0 

h* < h < h c Finite 

Finite 

0 

o 

A 

1 

4 

Finite (=£ -1/2) 

-1 


case of small a, we obtain L z /N ^ 0 for h < A. To see this, 
we evaluate the triplet pairing amplitudes 

(c^c^) = - ie~ l<p (u k nVku + Kmvjfcii), (10a) 

( c ki c ki) - _ ie l(f, (uk2iVk22 + Uk22Via\), (10b) 

whose pairing symmetries are (p + ip) wave. Figures [Ja-c) 
show the h dependence of the integrated pairing amplitudes 

f(p IP ) - J* 2 (UkllVkU + Mjfci 2 Vjfcii), (11a) 

F (p+ip) = J ^ 2^2 ( u k2\Vk22 + Uk22Vk2l), (Hb) 

for different a. For h < h*, both F are suppressed as 
seen in Fig.|4ja). 

On the other hand, for h* < h < h c , s wave is no longer 
dominant although it is topologically trivial. In Figs. Q a - C )> 
the ( p - ip)- wave pairing amplitudes E {p ~ ip ^ increase. Owing 
to the intraband (p + /p)-wave pairing potentials, the lower 
and upper bands acquire the orbital AM ±(€ z ) per electron, 
respectively. When the total numbers of electrons in the lower 
and upper bands are denoted by N+, respectively, the bulk or¬ 
bital AM is represented by L z = (€ Z )(N- - N+), and hence 
L z /N = (4 )(N- - N+)/(N- + N+) negatively increases as a 
function of h. Note that AL > N+ and (4) < 0. Especially, in 
the case of small a and A, we obtain L z /N = (£ z )h/p because 
of N + = m(p ± h)6(p ± h)l2n. 

For h> h c , since the upper band goes away from the chem¬ 
ical potential, we obtain N+ = 0 and L z /N = (€ z ). However, 
it is not correct to expect (4) = -1/2 just because the pairing 
potential is (p - ip) wave. In the presence of the Rashba SOI, 
not the orbital AM 4 = —id# but the total AM j z = l z + <t z /2 
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FIG. 2. (Color online) h dependence of the Chern number C (blue broken line), the total spin per electron SJN (red dashed line), and the bulk 
orbital AM per electron L z /N (black solid line). Black filled squares indicate L z /N obtained by numerical calculations in the circular disk. We 
set (a) a = 0.1, (b) a = 0.25, (c) a = 0.5, and (d) a = 1. 



a a 


FIG. 3. (Color online) a dependence of C, S z /N, and L z /N. The legends are the same as in Fig. [2] We set (a) h = 0.33 < A, (b) A < h = 0.5 < h c , 
(c) h = 0.67 > h c , and (d) h = 1 > h c . 


is conserved, i.e., [ j z ,Hj] = 0. We find that the wave func¬ 
tions Eqs. ( |2c| ) and ( |2d| are the spontaneous eigenstates of 
j z whose eigenvalues are ±1/2. Therefore, the bulk total AM 
J z = L z +S z always vanishes as seen in Figs.[2]and[3] For small 


a , spins are polarized in the z direction owing to the Zeeman 
interaction, which leads to (€ z ) - -1/2. Correspondingly, the 
(p + ip )-wave pairing amplitude F^ p+ip ^ is almost suppressed 
in Figs.[4ja,b). As a increases, spins are tilted, and (€ z ) goes 
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FIG. 4. (Color online) (a-c) h dependence of the pairing amplitudes F^ p ip) (blue broken line) and F^ p+ip) (red dashed line). We set (a) a - 0.1, 
(b) a = 0.5, and (c) a = 1. (d) h dependence of the normalized difference {F^ p ~ ip) - F^ p+ip) ) / {F^ p ~ ip) + F^ p+ip) ) for different a. 


to zero. Although E^ p+ip) itself is no longer suppressed in 
Fig. He), we find that the normalized difference between the 
pairing amplitudes {F^ p ~ ip) - F^ p+ip) ) / (F^ p ~ ip) + E^ p+ip) ) de¬ 
creases with a increasing as seen in Fig.[4jd). 


IV. DISCUSSION 

Here we discuss the effects of nonmagnetic impurities 
in time-reversal-broken topological SCs. In the topological 
phase, a strong nonmagnetic impurity behaves like a vortex 
core, which gives rise to the midgap bound statd 31 ! Since a 
vortex carries the orbital AM, such behavior seems plausible 
in the topological phase with the bulk orbital AM. The ratio 
\Eb/E g \ monotonically decreases as h increases 3 ^, where E b 
and E g represent the energy level of the bound state and the 
energy gap, respectively. This nonuniversal impurity effect 
can be characterized by the nonuniversal bulk orbital AM. 

To check this conjecture, we introduce a single impurity 
potential Vo 0(n - r) at the center of the circular disk and 
solve the BdG equations. Here 6(x) is the step function, and 
we set the radius of the impurity potential r x - 0.5 and the 
strength V 0 = 100. Remember that we set the number of the 
grid points N g = 1000 and the disk radius r c = 100. Fig¬ 
ures pla,b) show the h dependence of the eigenvalues E” =0 
for a = 1 and a dependence for h = 1, respectively. We 
can clearly see the midgap bound state in addition to the edge 
state E = 0 for h > h c . Since the impurity strength is suf¬ 
ficiently large, this bound state is what is called the “univer¬ 
sal midgap bound state” in Ref. [31 j Apparently, the obtained 
midgap bound state is not universal. The “universal midgap 
bound state” E ^ A 2 /£p in a topological SC 31 'relied on the 


edge state E(n) = -(n + 1/2)Ao/&f£ in a chiral /7-wave SC 22 . 
The former Ao is the s-wave pairing potential, but the latter is 
the chiral /7-wave one. We should replace Ao by A e ff, which 
depends on a and h. However, we cannot simply interpret A ef f 
as the (/7 - //7)-wave pairing potential because there coexists 
the (/7 + //7)-wave component as shown in Fig. [4] Although 
we assume the uniform pairing potential and do not solve the 
BdG equations self-consistently, we find that the energy level 
of this bound state \E^/E g \ and the bulk orbital AM 1 - \L Z /N\ 
are well correlated as seen in Figs [5jc,d). A self-consistent 
calculation may improve this correlation but is a future issue. 

With the appearance of the midgap bound states, the reduc¬ 
tion of the critical temperature T c also becomes relevant 32 . 
This crossover of the robustness against nonmagnetic impuri¬ 
ties is consistent with the crossover of the bulk orbital AM we 
have found. Moreover, we expect that the Rashba SOI makes 
the topological SC robust because it suppresses the bulk or¬ 
bital AM. 

We can also explain the vorticity-dependent impurity ef¬ 
fects on the vortex bound statJ^. As already mentioned in 
Sec. |TJ the bound state in an antiparallel vortex is more robust 
against nonmagnetic impurities than that in a parallel vortex. 
Furthermore, the low-energy scattering rate in the case of an 
antiparallel vortex is suppressed for small a as in chiral p- 
wave SCs, while it is similar to that in the case of a parallel 
vortex for large a as in s-wave SCs. The orbital AM carried by 
the antiparallel vortex is almost canceled by the intrinsic or¬ 
bital AM for small a, but not for large a. For sufficiently large 
a , the intrinsic orbital AM is totally suppressed, and two types 
of vortices cannot be distinguished. Therefore, we conclude 
that the impurity effects are characterized not by the integer 
Chern number but by the nonuniversal bulk orbital AM. 
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FIG. 5. (a) h and (b) a dependences of the eigenvalues E* 0 for the BdG equations with a single impurity, (c) h and (d) a dependences of the 
bulk orbital AM per electron 1 - \L Z /N\ (solid line) and the energy level of the midgap bound state \E b /E g \ (cross). We set (a,c) a = 1 and (b,d) 


h= 1. 


Although it may be difficult to measure the bulk orbital AM 
directly, there are two alternatives. One is the Hall viscosity, 
which is equal to half the bulk orbital AM in two-dimensional 
gapped systems at zero temperature 35 36 . The other is the spin 
magnetization, whose absolute value is equal to that of the 
bulk orbital AM but sign is opposite as discussed above. Es¬ 
pecially, in Fig.[2j the slope of S z /N at h —> 0 is nothing but 
the spin susceptibility^ = dS z /Ndh. In the presence of the 
SOI, Xzz is nonzero due to the interband Van Vleck spin sus¬ 
ceptibility. This is an indirect but simple way to measure the 
bulk orbital AM. 


V. SUMMARY 

To summarize, we have calculated the bulk orbital AM 
in the two-dimensional Rashba+Zeeman+v-wave topological 
SC. The Berry-phase formula coincides with numerical calcu¬ 
lations in the circular disk. Prior to the topological phase tran¬ 
sition at h = h c , we have found the crossover from v wave to p 


wave at h = h*, where L z /N starts to increase negatively. We 
have also found that the Rashba SOI tends to suppress S z /N 
and L z /N because the bulk total AM always vanishes. For 
large a, even in the topological phase, L z /N does not reach 
-1/2, which is the intrinsic value in chiral /7-wave SCs. We 
have explained the impurity effects on the topological SC: the 
appearance of the midgap bound state d 31 !I 33 1 the reduction of 
the critical temperature 32 , and the vorticity dependence of the 
vortex bound state 34 . 
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